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Abstract 

In this contribution to the conference "Beyond Einstein: Historical 
^H^i Perspectives on Geometry, Gravitation and Cosmology in the Twen- 

tieth Century" , we give a critical status report of attempts to explain 
the late accelerated expansion of the universe by modifications of gen- 
^ . eral relativity. Our brief review of such alternatives to the standard 

OO I cosmological model addresses mainly readers who have not pursued 

' the vast recent literature on this subject. 

in 

O ■ 1 Introduction 

OO 

^ ' The phenomenologically very successful cosmological 'concordance model', 

within the framework of general relativity (GR), leaves us with the mystery 
^ ■ of dark energy (DE). Since no satisfactory explanation of DE has emerged so 

it is certainly reasonable to investigate whether possible modifications 
of GR might change the late expansion rate of the universe. After all, GR 
has not yet been tested on cosmological scales. 

Modified gravity models have to be devised such that they pass the strin- 
gent Solar System tests, and are compatible with the rich body of cosmolog- 
ical data that support the concordance model (ACDM model). At the same 
time, the theories should be consistent on a fundamental level. Since we are 



*Invited talk at the conference "BEYOND EINSTEIN", Mainz, 22 - 26 September, 
2008. 

^See, e.g., [1], [2], and references therein. 
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dealing with higher spin equations, possible acausalities are, for instance, a 
serious issue. 

Apart from all that, one should not forget that the old profound vacuum 
energy problem [2] and the cosmic coincidence problem remain, and thus ex- 
treme fine tuning is unavoidable. This holds, of course, also for all dynamical 
models of DE [Ij. 

In my brief review I shall mainly concentrate on so called f{R) gravity. 
This is the simplest modification. Moreover, there have been some recent 
developments that I find interesting. After some generalities, many of you 
know very well, I shall discuss the weak field limit and Solar System tests. 
For some time there was a lot of confusion on this issue, with conflicting 
statements, but was eventually clarified. We shall, however, see that the weak 
field approximation may break down, and a so-called Chameleon mechanism 
can be at work that hides a scalar degree of freedom of the theory on solar 
system scales. 

There are f{R) models that pass the solar system tests and are cosmolog- 
ically almost indistinguishable from the successful ACDM model. Recently 
it was, however, discovered by Kobayashi and Maeda that these models are 
in serious trouble in the strong-field regime. 

Some of the other modified gravity theories are even in greater difficulties. 
This will be briefly discussed in a final part. 



The simplest possibility of modifying GR is to replace the Einstein-Hilbert 
action i? — 2A of gravity by a nonlinear function f{R) of the Ricci scalar 
This introduces an additional scalar degree of freedom that can lead to 
an accelerated expansion of the universe at late times, induced by the Ricci 
scalar. One may call this "curvature DE" or "dark gravity" . The function / 
is quite arbitrary and the theory looses, of course, a lot of predictive power. 
As many other people, I regard this class of modified gravity theories as 
instructive phenomenological toy models, that change gravity in the infrared. 
The variation of the gravitational action is: 



Diffeomorphism invariance implies that the tensor in the curly bracket has a 
^For an extensive review and literature, we refer to [3]. 



2 Metric f(R) gravity 




V^Vpf{R)}5g''^V^ d'x . (1) 
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vanishing covariant divergence. Therefore, the field equation 



+ go.pUf\R)-V^Vpf\R) = STiGT^p (2) 

imphes that the energy-momentum tensor T"^ is divergence-free. This im- 
phes, by a general result of Hawking, that matter propagates causally if T"^ 
satisfies the dominant energy condition. As expected from Lovelock's the- 
orem, the field equation is of fourth order in the metric if / is not a linear 
function. 

It is easy to see that the de Sitter or anti-de Sitter metric, with = 
Ag^is, is a vacuum solution, if there is a constant A satisfying /(4A) = 
2Af (4A)| This indicates that the theory may naturally lead to cosmological 
acceleration. More on this later. 

At first sight one may think that experience and insight from GR may not 
help us to get some understanding of what the complicated fourth order field 
equations may describe. It is, however, known since long that f{R) gravity 
models can be reformulated as scalar-tensor theories [4j, [S]. To show this, 
we first rewrite the field equations in the following form (k^ := SttG) 

f'{R)G^f3 = K^T^p + ]^g^p[f{R)-Rf\R)] 

- gap^f\R) + VoSpf\R) ■ (3) 

It is natural to introduce the scalar field := f'{R)- We assume that f" ^ 0, 
so that /' is at least locally invertible: f oTZ = id. Furthermore, let the U 
Legendre transform of /: 

w(0) = 7^(0)0-/(7^(0)) (4) 

(thus W = R). With this we can rewrite (3) as 

(j)GaP = n'^TafS - ^ga/3l^{(p) - [gaf304> " V^V/J^]. (5) 

This is just the Brans-Dicke equation with the Brans-Dicky parameter ubd = 
plus a potential termQ This indicates that the weak field limit may be in 



■^If this transcendental equation has a solution, any vacuum solution of GR with the 
corresponding A is obviously a vacuum solution of (2). 

^Therefore, one expects that the Cauchy problem is well-posed. This is certainly the 
case for the vacuum theory, but with matter the problem is not completely settled (see 
[5], [7]). In [7] two first order strongly hyperbolic formulations of scalar-tensor theories 
are presented, which however do not include the exceptional case uj = —3/2. 
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conflict with solar system tests, because these imply that the parameter ubd 
has to be very large: ubd > 40'000. We shall see that this is indeed the case, 
but thanks to the potential term there is an interesting way out. 
Taking the trace of the original field equation (3) we obtain 



3nf\R) + Rf\R)-2f{R) = K'T. 
In terms of the scalar field d> this becomes 



3n0 + 2W(0) 



K^T. 



(6) 



(7) 



This nonlinear scalar field equation will play a crucial role. It shows that the 
scalar degree of freedom is truly dynamical. In contrast to GR, the scalar 
Ricci curvature does no more track the matter distribution. 

We note at this point, that the field equations (5) follow from the following 
action 

S = ^j[<PR- d'x + Sm. (8) 

Since no kinetic energy for the (p-field appears in this action, one may be 
tempted to conclude that (p is not a dynamical field, but we have seen that 
this is not the case. 

For certain problems it can be useful to pass to a mathematically equiv- 
alent description by performing the conformal transformation (first studied 
by Pauli in letters to Jordan in 1953 [H]): 



g^u = exp 



exp 



KLf 



(9) 



In terms of the new metric and the scalar field ip, called the Einstein frame, 
the action becomes 



S 



EF 



where (3 



and 



g d x + Suig 



(10) 



2k^ 



-2/3^ [e^^7^(e^^) - / (7^(e^^))] . (11) 



In contrast to the original Jordan frame description, the gravitational part of 
the action takes a canonical form familiar from GR, but the coupling to mat- 
ter is non-minimal. This implies that relative to the Levi-Civita connection 
belonging to the metric g^iy, the energy-stress tensor is no more conserved. 
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In the Einstein frame matter feels a new 'fifth force' due to gradients of (p. 
While Newton's constant is everywhere the same, the local particle physics 
thus varies. In the Jordan frame, on the other hand, the laws of physics in 
local inertial frames are universal, but the effective gravitational "constant" 
(G/0) becomes space-time dependent. Since the two descriptions are math- 
ematically equivalent, observables are frame independent. It is then just a 
matter of convenience which one prefers to use. In what follows we will 
always work in the Jordan frame, except at one instance. 

3 Generalized Friedmann models 

It is straightforward to derive the modified Friedmann equations. We con- 
sider only Friedmann-Lemaitre (-Robertson- Walker) spacetimes with are spa- 
tially fiat. If a{t) denotes the scale factor and H = a/a the Hubble rate, one 
finds 



are effective fluid contributions due to curvature ('curvature dark energy'). 
The sign of the corresponding effective equation of state parameter Weff ■ = 
Peff/peff is determined by that of Pe// since pe// has to be non- negative. 
Simple choices for f{R) give strongly negative values for w^ff. For example, 
f{R) = R — / R gives Weff = —2/3. As we noted before, the effective 
gravitational "constant" in (12) and (13) is G/f, and is thus i?-dependent. 

In general none of the standard energy conditions is satisfled for f{R) 
models. In particular, | Pe// | does not have to be smaller than p. 

Clearly, the 'energy conservation' p = — 3(p + P) follows, as in GR, from 
the fleld equations, i.e., from (12) and (13). 

Since a/ a = H + we may regard the evolution equations (12) and 
(13) as a system of equations for H and R. This dynamical system has 

^For symmetry reasons T^" has the form of an ideal fluid 




(12) 



(13) 



Peff = l{Rf'-f)-3HRf", 

P,/; = R^f" + 2HRf" + Rf" + ]^{f-Rf) 



(14) 



(15) 
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been extensively studied (see, e.g., [^). It is a mathematical fac1@ that 
for any given expansion history a{t) there exists a (non-unique) function / 
that reproduces this history by the corresponding f{R) model. This does, 
however, not guarantee that the sequence of radiation-matter-acceleration 
eras is also reproduced. Indeed, the analysis in [10| shows that non-linear 
f{R) models that reproduce, for example, exactly the history a(t) of the 
ACDM model, do not have the right sequence of cosmological eras with 
the required density parameters. However, it is possible to reproduce this 
sequence of eras if some deviations from the given history H{z) are tolerated. 
An example for this, that is compatible with current observations, has been 
given by Hu and Sawicky pTj . 

Since this model will also later play a role, we present it here. Its analytic 
form reads 

f(R] = R — m — -— - — -T , m := 16 

(see Fig. 1), with suitably chosen parameters ci,C2 {n is a positive integer, 
and po denotes the present average cosmic density). The corresponding his- 
tory leads to a curvature equation of state parameter Weff{z) that deviates 
from -1, but these deviations can be kept sufficiently small. At high redshifts 
Weff{z) becomes smaller than -1, a possibility that can be checked with fu- 
ture observations. Such a crossing of the so-called phantom line would be 
interesting. 

We do not discuss here the evolution of linear cosmological perturbations 
for f{R) models, that may imply interesting testable deviations from the 
concordance model. (For some references on this, see the last section.) 



4 Weak field limit for spherically symmetric 
sources 

We remarked earlier that for f{R) gravity the Schwarzschild or Schwarzschild- 
de Sitter metric is often a vacuum solution, e.g., for the model (16). This 
does, however, not guarantee that the theory passes solar system tests. This 
would be the case if this vacuum solution can be matched to an interior solu- 
tion (as Schwarzschild showed for GR). We shall see that this is generically 
not possible. In this section we address this issue in the weak field limit. 

As a preparation we linearize the scalar field equation (6) about a back- 
ground de Sitter universe with A = -Rs/4, where the Ricci scalar Rg satisfies 

f'.Rs -2fs = (17) 

^For a dynamical system analysis of this reconstruction, see [lOj . 
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Figure 1: Plot of | {f{R) — R)/m? \ from Hu and Sawicki [TT]. (In the graph 
/(i?) denotes our f{R) — R, and Jrq is our f'{R) — 1 today.) The ratio Ci/c2 
is fixed such that f{R) approaches for large R/m? the ACDM value. The 
value of R/m? at the present epoch is about 40. Thus only the part of the 
graph to the right of this value is relevant. 
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ifs '■= f{Rs), etc). Let R{r) = Rg + SR{r), and linearize the trace equation 
for a local source: 

3f:n6R + {RJ'J - Q5R = k^T. (18) 

Since T vanishes for the background, we regard it as of first order. The last 
equation shows that the scalar field 5R has an effective mass given bjlll 

2 _ f's " f's'Rs _ (/s)^ ~ 

After considerable confusion, it was shown in [T3j that the Eddington- Robertson 
parameter 7 is not equal to 1, as in GR and also observationally to high ac- 
curacy, but 7 = 1/2, if the following conditions are satisfied: 

(i) Linearization of f{R), f\R) about Rg is allowed. 

(ii) nRs) ^ 0. 

(iii) The Compton wave length l/m^ is much larger than the size of the 
solar system. 

(iv) The deviations of the gravitational field from the de Sitter background 
metric can be treated in first order. 

Remarks. These conditions are not always satisfied. If f"{Rs) = 0, then 
7 = 1 as in GR. Condition (iii) can be violated, for instance by fine tuning 
the parameters in 

fiR)=R+^R^-fiyR. 

The only way to escape the destructive consequence 7 = 1/2 and maintain 
the late cosmic acceleration, is to invoke a "chameleon mechanism" for the 
scalar degree of freedom. 



5 Chameleon mechanism 

The chameleon effect was discovered by Khoury and Weltman [T3] in scalar 
field models of DE. Scalar fields with self-interactions may directly couple 
to matter as strong (or even stronger) as gravity and still satisfy all current 
constraints. The reason for this is that the effective mass of the scalar field 
depends on the local density. So there is the possibility that the Compton 
wave length is sufficiently small on Earth to satisfy all Laboratory bounds, 

'''It turns out [H] that the non-negativity of the expression in (19) is the stabihty 
condition of the de Sitter spacetime with respect to small inhomogeneous perturbations 
of the f{R) model (without matter). 
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while it is much larger in the Solar system and still much larger on cosmo- 
logical scales. 

For illustration, consider a scalar field model satisfying the non-linear 
equation 

□y, = K'/zlv'), Veff = V{v)-B{P^/Mp{)T, (20) 

where T is the trace of the matter part of the energy-momentum tensor {k. — p 
if the pressure can be neglected). The dependence of Veff on p can imply 
that d'^Vf.ff/dip^ at the effective minimum is much smaller for a low density 
background than in a high density environment (see Fig. 2). This density 
dependence can lead to a thin-shell effect: ip varies for a macroscopic body 
only over a thin surface layer, leading to a weak fifth force. This behavior is 
intimately related to the non-linear nature of Chameleon field theories. An 
equation of the type (20) is obtained for f{R) models in the Einstein frame. 
By transforming Eq. (7) one finds 

dip 2 

Only non-relativistic matter contributes to T. If we define p by 

T^-p=: -e'^^PI^^^p, 
then p is conserved in the Einstein frame. In terms of this quantity we obtain 

□V=^, f/e//(y^) = f/(y^)+e-(^/^)-p, 

which is of the form (20) with an exponential function B (as in Fig. 2). In 
what follows, we do not make use of this Einstein frame formulation. 

The possibility of a Chameleon mechanism for the scalar degree of freedom 
of f[R) gravity models has been studied in several papers, e.g., in [16], [T7] . 
We discuss here briefiy part of the work by Hu and Sawicky [llj^ . 

These authors study in the framework of f{R) models non-relativistic 
stars like the Sun with weak gravitational fields, but do no more linearize 
the equation for = f\R)- For a static situation Eq. (6) reduces to 

3V^f'{R) + Rf\R) - 2f{R) = K^T ^ -K^p. (21) 

It is a good approximation to replace the Laplacian of the slightly curved 
spatial metric by the fiat space Laplacian (for which we use the same symbol). 
Given a density profile p{r) from a solar model, Eq. (21) becomes a non- 
linear field equation for R (or (fi). Hu and Sawicki choose the model (16) and 



^For a simplified discussion in the Einstein frame, see 
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pexp(|3(|)/ M ) 



V((^) 



Figure 2: Typical effective potential of the form (20), whose density depen- 
dence can lead to a Chameleon effect (from ref. |15j). 



impose the following boundary conditions: Deep inside the star /' assumes 
the value with R = K^p (implied by GR). Very far away (r = IO^Vq) the outer 
boundary condition /' = f'{R = n'^Pg) is imposed, where Pg is the average 
galactic density in the solar vicinity [pg = lO"^"^ g cm~^). These boundary 
conditions correspond approximately to the minima of the effective potential 
belonging to (21) At this point we consider them as part of the model. 
The chosen density profile is shown in Fig. 3 (solid line), while the numerical 
solution of the boundary value problem for R{r) is shown by the dashed line. 
(The parameter fjiQ in this figure is /' — 1 for the present average cosmic 
scalar Ricci curvature.) A blown up version of the region where R does not 
track the GR limit /t^p outside about 1 AU is shown in Fig. 4. 

Once p(r) and R{r) are known, the field equations (3) determine the 
Einstein tensor, from which the metric in the weak-field limit can easily be 
computed (Poisson integrals). 

Since R deviates from the GR value /t^p only in a very low density shell, 
the Solar System tests present no problem. For example, the 7 parameter 
becomes 

^ ~ 1 2U — (22) 



^The effective potential is defined by = -^p+\{2f - Rf) 
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Figure 3: Realistic density profile of the solar interior and vicinity (solid 
curve), and solution R/k^ of the scalar field equation for n = 4 and field 
amplitude |/'(i?o) ~ 1| = 0.1, with boundary conditions decribed in the text 
(dashed line). (From with the same change of notation as in Fig. 1.) 
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Figure 4: Blown up version of Fig. 3 for various choices of the cosmological 
field amplitude f'{Ro) — 1 (equal to Jrq in the notation of [llj). 
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where M is the total mass of the star, and 

Me// = 47r j {p- R/n^y^dr. (23) 

For the solutions shown in Figs. 3,4 Me// <C M, thus 7 1 to high accuracy. 
So this looks very good, but it has to be emphasized, that this conclusion 
rests on the assumption, that the galaxy is in the high- curvature regime [R ~ 
K^p) with respect to its own density profile. The validity of this assumption 
depends, as Hu and Sawicki stress, "on both the structure of the galactic 
halo and its evolution during the acceleration epoch". This dependence is 
irritating, but reflects the fact that there is no Birkhoff theorem for f{R) 
gravity models. To decide about the proper boundary conditions one would 
have to study - horribile dictu - galaxy formation for f{R) models in N-body 
simulations. 

We shall see in the next section that the Chameleon mechanism can only 
work if the star has a surrounding medium, e.g., dark matter. 



6 Nonexistence of relativistic stars in f{R) 
gravity ? 

Recently it was discovered by Kobayashi and Maeda [l9j that models which 
incorporate the Chameleon mechanism cannot describe relativistic stars. 
This important result presumably excludes f{R) gravity as a viable mod- 
ifications of OR. In this section we briefly describe the content of the paper 
by Kobayashi and Maeda (abbreviated as KM). 

Especially for the numerical part of their work KM use the following 
model 

f{R) = R + XR'[{l + ^)--l\ {X,R>0), (24) 

that was adopted by Starobinsky [20] to incorporate the Chameleon mech- 
anism. Since only the qualitative behavior of the potential V{(p) := U/cf^ 
near the de Sitter minimum matters, the results will also apply to (16) and 
other models [21]. We set the de Sitter value Rg = Xsm, then A is uniquely 
determined by Xs (for a given integer n). The same holds for the de Sitter 
value (ps of (j). We note that 

0(i?) = l-2nA^(l + |^) , (25) 



^"Note that V{(j)) is closely related to U in (11). It is easy to see that the de Sitter value 
for R is mapped to the value of 0, where V takes its minimum. 
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which shows that a curvature singularity {R — > ±00) is mapped to = 1. 

It should be remarked that V{(f)) is a multivalued function, but this is no 
worry because only a relatively small interval about 0s really matters. 

For this model KM study spherically symmetric stars. The metric is 
parameterized in Schwarzschild coordinates as 

g = -N(r)df + -^dr^ + r\de^ + r^sin^ 9 dip^). (26) 
B{r) 

It is easy to generalize the GR structure equations to nonhnear f{R) models. 
The (tt) and (rr) components of the field equations (5) become 







^{-1+B + rB) = -SttGI^-^V-S 0,,+ + ,(27) 



T 



' 2 Bt 



^A-^ + B + rB^) = 8.GP-<f>W-B{l + §^)<p.. (28) 

An index r denotes differentiation with respect to r. The scalar field equation 
(7) becomes 



B 



,2 Nr Br 



SttG . ^ 2(f)' 
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+ (^ + 2i + ii)'^^ =^(-p + 3P) + ^Vm (29) 



3 ' ' '3 



Recall that the energy-momentum 'conservation' gives, as in GR, 

N 

Pr+^{P + P)^0- (30) 

We turn to the boundary conditions. Near the center of the star at r = 
the various functions are expanded in powers of r^, making also use of the 
scaling freedom of the time coordinate. Some of the coefficients are deter- 
mined in terms of others by the field equations. So far KM have assumed, for 
simplicity, that the energy density of the star is constant, thus generalizing 
Schwarzschild's interior solution. Given p and the central values Pc, (pc of P 
and 0, the basic Eqs. (27)- (29) can be integrated outwards from the center 
to the surface r = i?* of the star, which is defined by P{R*) — 0. Note that 
Eq. (30) gives 

r + P 1 2 

N(r)^ ^ % . (31) 

Prom the surface the vacuum equations are integrate to sufficiently large r. 
The starting values at the center are then varied until the de Sitter solution 
with ^ as r ^ cxD is assumed. 

It turns out that there are no solutions if the gravitational potential $ := 
{l—N)/2 is larger than some value ^maxi which is typically about 0.1. (Recall 
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that for GR ^max = 4/9.) When KM tried to find solutions with larger 
the Ricci scalar diverged. Thus, for the studied class of f{R) models, there 
are no stars with strong gravitational fields. The authors give also 
analytic arguments for this conclusion, that was originally suggested by [22] . 
These are based on a mechanical interpretation of the scalar field equation 
(29). 

For non-relativistic stars it turns out that the thin-shell condition is vio- 
lated, when there is no matter outside the stellar surface, and therefore the 
parameter 7 is close to 0.5. For such stars one can also derive good analytic 
approximat ions . 

In the previous section we saw that the Chameleon effect can work for non- 
relativistic stars if surrounding matter is taken into account. Surrounding 
matter does, however, not change the non-existence statement for strong 
gravitational fields, as is shown by KM. 

In view of these results f{R) gravity models are in serious trouble. More 
realistic equations of state will presumably not change this conclusion. 

7 Inclusion of other curvature invariants 

There are a number of studies [23], that include other curvature invariants, 
such as R^^R'^^, Rap^sR"^^^ . Such models are in most cases unstable, like 
mechanical Lagrangian systems with higher derivatives [24j"1. An exception 
seem to be Lagrangians which are functions of R and the Gauss-Bonnet in- 
variant Rgb = R^ — 4,R^uR^^ + Raji^sR"^^^ ■ By introducing two scalar fields 
such models can be written as an Einstein-Hilbert term plus a particular 
extra piece, containing a linear coupling to Rgb- Because the Gauss-Bonnet 
invariant is a total divergence the corresponding field equations are of sec- 
ond order. This does, however, not guarantee that the theory is ghost-free. 
In [26] this question was studied for a class of models [23] for which there 
exist accelerating late-time power-law attractors. It turned out that in a 
Friedmann background there are no ghosts, but there is instead superluminal 
propagation for a wide range of parameter space. This acausality is remi- 
niscent of the Velo-Zwanziger phenomenon [27] for higher (> 1) spin fields 
coupled to external fields. It may very well be that it can only be avoided 
if very special conditions are satisfied. Ghosts of Gauss-Bonnet cosmologies 
have also been studied in [28] . 

In addition to these problems, it appears unlikely that the devastating 
difficulties we have encountered for f{R) models will disappear when other 

^^This paper contains a discussion of a generic instability of Lagrangian systems in 
mechanics with higher derivatives, that was discovered by M. Ostrogradski in 1850 [5S]. 
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curvature invariants are included. 



8 First-order (affine) modifications of GR 

The disadvantage of complicated fourth order equations can be avoided by 
using the Palatini variational principle, in which the metric and the sym- 
metric affine connection (the Christoffel symbols T"^,^) are considered to be 
independent fields!^ 

For GR the 'Palatini formulation' is equivalent to the Einstein-Hilbert 
variational principle, because the variational equation with respect to T°'^^ 
implies that the affine connection has to be the Levi-Civita connection. 
Things are no more that simple for f{R) models: 



S 



2k 



f{R) + L 



matter 



(32) 



where R[g,V] = g^'^RajslT], Ra/3[T] being the Ricci tensor of the independent 
torsionless connection T. The equations of motion are in obvious notation 



f'{R)R(,,)[r]--f{R)g,, 
Vl [V^gfW) = 



(33) 
(34) 



For the second of these equations one has to assume that Lmatter is function- 
ally independent of F. (It may, however, contain metric covariant deriva- 
tives.) 

Eq. (IMl) implies that 



V 



99 







(35) 



f{R)9tiu- Hence, the F"^^, are 



for the conformally equivalent metric g^i, 
equal to the Christoffel symbols for the metric g^j^iy. 
The trace of (l33|) gives 

Rf'iR) - 2f{R) = kT. 

Thanks to this algebraic equation we may regard i? as a function of T : R = 
TZ{T). In the matter-free case it is identically satisfied if f{R) is proportional 
to R^. In all other cases R is equal to a constant c (which is in general not 



^^This approach was actually first introduced by Einstein 29J. This is correctly stated 
in Pauli's classical text on relativity (p. 215). 
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unique). If /'(c) 7^ 0, eq. flM|) implies that T is the Levi-Civita connection 
of g^^, and flH^ reduces to Einstein's vacuum equation with a cosmological 
constant. In general, one can rewrite the field equations in the form of 
Einstein gravity with nonstandard matter coupling^. 

- ^[V;./V./'-^^^.(V/')1 + (V^V.-^?^.n)/'. (36) 

With this reformulations it is, for instance, straightforward to develop cos- 
mological perturbation theory [31]. 

For some time this modification of GR looked promising. But now one 
can ignore it because of the following major drawbacks. 

1. Since the vacuum theory is identical with GR including a cosmological 
constant, the metric g outside a spherically symmetric star has to be the 
Schwarzschild-de Sitter metric. This does, however, not guarantee that the 
Solar System tests are satisfied. For this we have to know whether there 
are interior solutions that match the exterior metric field. This was studied 
in several papers, with a negative result for physically relevant equations 
of state. Technically it was shown in [32] that for polytropic equations of 
state with an adiabatic index F in the interval 3/2 < F < 2 true curvature 
singularities develop. From this one can guess with confidence, that the 
theory cannot describe white dwarfs - for example. 

2. For nonlinear Palatini f{R) gravity the Cauchy problem is not well- 
posed [33pl . 

Both of these unacceptable shortcomings have in the final analysis a com- 
mon origin. The good thing about the field equations in the form (36) is that 
it they are of second order in the derivatives of the metric. What leads to the 
mentioned difficulties is that the right hand side of (36) is at least of second 
order in the matter variables, because i? is a function of T. Apart from ideal 
fluids, T is usually of flrst order, and then the right hand side is even of third 
order in the matter flelds. This induces locally sensitive dependence of the 
metric on rapidly changing matter flelds, in contrast to GR (and Newtonian 
gravity) where such dependencies are smoothed out. 

"'^^It is shown in |30j that if the matter action is independent of F, the theory is dynam- 
ically equivalent to a Brans-Dicke theory with the special Brans-Dicke parameter —3/2, 
plus a potential term. 

^^In ,33J it is shown that the basic system of equations in vacuum can not be rewritten 
as a system of only first order, since can not be eliminated, except of course if \3(f) = 
(e.g., for the vacuum theory). 
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9 Concluding remarks 



A positive aspect of the largely negative outcome of the previous discussion 
seems to me that the distinguished role of GR among large classes of gravity 
theories has once more become apparent. We know, of course, that GR 
is an effective theory, and that quantum theory will produce all sorts of 
induced terms (a phenomenon that is well-known from QED), but stopping 
any expansion after a few terms will hardly lead to a consistent theory that 
agrees with observations on all scales. 

Some of the modified gravity theories, such as f{R) or braneworld models, 
may perhaps be of limited use for testing GR on cosmological scales. Guided 
by such model^, there have recently been some interesting attempts to de- 
velop a parameterized post-Friedmann description of gravity that parallels 
the parameterized post-Newtonian description of Solar System tests (dis- 
cussed by C. Will at this meeting). In contrast to the latter, there appear 
unavoidably some free functions, instead of just a bunch of parameters, in 
the description of the evolution of inhomogeneities [37], [3H]. It will, there- 
fore, be difficult to discriminate between dark energy and modified gravity, 
but this remains a major goal for years to come. One can hope that this will 
eventually become possible with better data on the CMB background, weak 
gravitational lensing, and the growth of large scale structures. 
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